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TREND-FREE  BLOCK  DESIGNS:  THEORY1 
by 

Ralph  A.  Bradley  and  Ching-Ming  Yeh 
Florida  State  University 

A common  polynomial  trend  in  one  or  more  dimensions  is  assumed  to  exist 
over  the  plots  in  each  block  of  a classical  experimental  design.  An  analysis 
of  covariance  model  is  assumed  with  trend  components  represented  through  use  of 
orthogonal  polynomials.  The  objective  is  to  construct  new  designs  through  the 
assignment  of  treatments  to  plots  within  blocks  in  such  a way  that  sums  of 
squares  for  treatments  and  blocks  are  calculated  as  though  there  were  no  trend 
and  sums  of  squares  for  trend  components  and  error  are  calculated  easily.  Such 
designs  are  called  trend-free  and  a necessary  and  sufficient  condition  for  a 
trend-free  design  is  developed.  It  is  shown  that  these  designs  satisfy  optimality 
criteria  among  the  class  of  connected  designs  with  the  same  incidence  matrix. 

The  analysis  of  variance  for  trend-free  designs  is  developed.  The  paper  concludes 
with  two  examples  of  trend-free  designs. 
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1 . Motivation,  introduction  and  summary.  Various  experimental  designs  have 
been  considered  for  the  two-way  elimination  of  heterogeneity,  the  latin  square 
being  the  best  knovm.  Many  experimental  situations  arise  in  which  response  may 
be  affected  by  the  spatial  or  temporal  position  of  the  experimental  unit  or  plot 
within  a block  and,  frequently,  an  assumption  of  a common  polynomial  trend  of 
specified  degree  over  plots  within  blocks  may  be  appropriate.  The  v * v latin 
square  may  be  used  to  eliminate  the  effects  of  a common,  one-dimensional, 
polynomial  trend  to  degree  v - 1 over  the  nlots  in  rows  of  the  square.  Ques- 
tions of  criteria  for  and  the  optimality  of  other  block  designs,  complete  or 
incomplete,  that  eliminate  the  effects  of  common  polynomial  trends  of  specified 
degree  over  the  plots  within  blocks  are  investigated.  Examination  of  the  exist- 
ence and  construction  of  the  desired  designs  is  deferred  to  a subsequent  paper. 

Experimental  designs  to  be  used  in  the  presence  of  trends  to  avoid  the 
complications  of analysis  of  covariance  and  to  increase  design  efficiencies  have 
been  developed  by  others.  Cox  (1951,  1952,  1958)  considered  the  assignment  of 
treatments  to  plots  ordered  in  space  or  time  without  blocking  and  with  a trend 
extending  over  the  entire  sequence  of  plots.  Box  (1952)  and  Box  and  Hay  (1953) 
in  similar  experimental  sequences  investigated  choices  of  levels  of  quantitative 
factors.  Hill  (1960)  combined  the  designs  of  Cox  and  Box  to  form  new  designs 
to  study  the  effects  of  both  qualitative  and  quantitative  factors  in  the  presence 
of  trends.  Daniel  and  V.'ilcoxon  (1966)  and  Daniel  (1976,  Chap.  15)  provided 
methods  of  sequencing  the  assignments  of  factorial  treatment  combinations  to 
experimental  units  to  achieve  better  estimation  of  specified  factorial  effects, 
again  in  the  presence  of  a trend  in  time  or  distance.  Phillips  (1964,  1968a,  1968b) 
illustrated  the  use  of  magic  figures  (magic  squares,  magic  rectangles,  etc)  for 
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the  elimination  of  trend  effects  in  certain  classes  of  one-way,  factorial,  latin- 
square,  and  gTaeco-latin-square  designs;  a single  trend  was  assumed  again  to 
affect  a'l  of  the  observations. 

We  consider  block  designs  to  compare  v treatments  in  b blocks  of  equal  size 
k sv  such  that  each  treatment  occurs  at  most  once  in  each  block.  A common 


polynomial  trend  is  assumed  to  exist  over  the  plots  in  each  block.  The  trend  may 
be  in  one  or  several  dimensions  and  is  expressed  in  terms  of  p s (k  - 1)  ortho- 
gonal components.  The  appropriate  classical  fixed -effects  model  for  the 
general  block  design  and  its  assumptions  are  used  with  the  addition  of  orthogonal 
terms  representing  the  assumed  trend.  The  problem  is  to  assign  treatments  to 
plots  within  blocks  so  that  the  known  properties  of  the  ordinary  analysis  of 
variance  for  treatment  and  block  sums  of  squares  are  preserved  and  variation 
due  to  the  trend  may  be  removed  from  the  error  sum  of  squares.  Mien  the  desired 
designs  exist,  we  shall  call  them  trend-free  designs.  We  shall  abbreviate  and. 


for  example,  use  TF^CB  for  a complete  block  design  free  of  a common,  one- 
dimensional trend  of  degree  p within  blocks  and  TF  BIB  and  TF  PBIB  for  similar 

p p 

balanced  and  partially  balanced  incomplete  block  designs.  Additional  subscripts 


may  be  added  if  the  trend  is  in  several  dimensions. 

In  this  paper,  necessary  and  sufficient  conditions  for  the  existence  of 
trend-free  block  (TFB)  designs  are  obtained,  necessary  analysis  of  variance  is 
developed, and  optimality  properties  are  demonstrated.  Several  examples  of  TFP. 


designs  are  included  but  results  on  design  construction  are  not  included. 


2.  Notation,  definitions  and  model.  Let  v treatments  be  applied  to  plots 
arranged  in  b blocks,  each  of  size  k,  k £ v.  Each  plot  receives  only  one 
treatment  and  each  treatment  occurs  at  most  once  in  a block.  The  plot  positions 
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in  a block  are  indexed  by  m-dimensional  vectors  of  positive  integers, 

t = (tj,  t^),  tu  * 1,  . ..,  su>  u = 1,  m,  where  su  is  the  number  of 

plot  positions  in  the  u-th  dimension,  TIs  = k.  The  polynomial  trend  extending 

u 

over  plots  in  each  block  is  a function  of  the  plot  position  t and  is  expressed 

as  a linear  function  of  m-dimensionai  orthogonal  polynomials  of  the  form, 

m 

<f>a(t)  = a .(t.,  ....  tm)  = n d>a  (t  ),  where  * (t  ) is  a one-dime:’ sional 

orthogonal  polynomial  of  degree  au  satisfying  the  orthogonality  conditions,  for 

U—  1,  . « . , m. 


(2.1) 


and 


(2.2) 


nU 


t =1 
u 


t>  (t  ) = 
a u 
u 


0 


L $ (t  ) ♦ , (t  ) - 
.L  . a u a 1 v u' 
t =1  u u 

u 


1 if  a = a' , 

u u 

0 if  a * a' . 

u u 


Note  that  multidimensional  orthogonal  polynomials  may  be  obtained  as  indicated 

m 

from  one-dimensional  ones  and  that  k = ns.  Tables  of  one-dimensional  ortho- 

u=l  u 

gonal  polynomials  for  equally  spaced  variables  are  given  by  Fisher  and  Yates  (1957). 

The  mathematical  model  is  a simple  extension  of  the  classical  model  for 
general  block  designs,  trend  terms  added,  and  is  written 

v • 

(2.3)  y..  = y ♦ 7 6..  t.  + 8.  ♦ 7 6 <J>  (t)  ♦ e.^, 

3%  jt  i 3 a V~  jt’ 

j - 1.  ....  b,  t - (t,,  ....  tj,  t„  - 1,  ...,  s.,  u = 1,  ...,  m. 


u' 


where  y^  is  the  observation  on  plot  position  £ of  block  j,  y,  and  6^  are 
respectively  the  usual  mean,  treatment  and  block  parameters,  A is  an  index  set 


4 


of  p,  m-dimensional , non-zero  vectors  of  the  form  a,  and  0 is  the  regression 

»v 

coefficient  of  <J>  (£) . The  trend  effect  on  plot  position  t is  £ 9 <j>  (t)  , not 


aeA 


s a 


dependent  on  the  particular  block  j.  Designation  of  the  treatment  applied  to 
plot  (j,  t)  is  effected  through  the  indicator  function, 

1 if  treatment  i on  plot  (j , t) , 

0 otherwise. 


(2.4) 


The  errors  e.  are  assumed  to  be  i.i.d.  with  zero  means  and  they  will  be  taken 
J £ 

2 

to  be  normal,  N(0,  a ) , in  discussions  of  distribution  theory  below.  The  model 
(2.3)  will  be  regarded  as  a fixed-effects  model  unless  specifically  stated  to 
be  otherwise  in  particular  remarks. 

Choices  of  values  of  specify  particular  block  designs.  The  construction 
of  TFB  designs  is  the  determination  of  values  of  the  6*  to  meet  criteria  to  be 
developed. 

Some  matrix  notation  will  be  needed.  Let  I be  the  identity  matrix  of 

~n 

order  n,  0 _ be  the  m x n null  matrix,  J „ be  the  m * n matrix  with  unit 
~m*n  'Tnxn 

elements,  1^  be  the  n-dimensional  column  vector  with  unit  elements,  and  BBC, 
the  Kronecker  product  of  matrices  B and  C.  Lexicographic  order  is  defined: 


Definition  2.1.  Two  distinct  vectors  of  non-negative  integer  elements, 
i * (ij,  ....  im)  and  i = ( j ^ , ...,  j^)  are  in  lexicographic  order  if  i is 
ordered  before  j,  whenever  ij  * j^,  ....  is  j = js_j»  is  < j for  some  s, 

1 £ s £ n.  Several  distinct  m-element  vectors  are  in  lexicographic  order  if 
all  pairs  of  vectors  are  in  lexicographic  order. 

The  model  (2.3)  in  matrix  notation  is 


(2.S) 


x ■ x 1 ♦ i - y * M * M * h6~  * £ 
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Y and  E being  bk-element  column  vectors  with  elements  y.  and  e.  in  subscript 

J k,  3 £ 

lexicographic  order  on  ( j , t ) , x'  = (x]L,  . ..,  xv),  £'  = (Bj,  ....  Bfe) , 

»’  - <ea  » •••»  ea  ).  X'  • (w.  I'.  &’).  X - XT,  x6,  Xg),  ^ = lbk, 

~1  ~p 

X^  = •••*  Xg  = Ib  ° ~k’  ~e  = ib  9 where  ~j  is  the  k X v matrix  with 

as  the  i-th  element  in  row  t and  £ is  the  k x p matrix  with  ^(t)  in  row  t 
and  column  a,  a e A,  and  the  rows  of  A^  and  £ being  in  subscript  lexicographic 
order  on  t.  For  convenience,  the  vector-elements  of  a are  taken  also  to  be  in 
lexicographic  order.  Conditions  (2.1),  (2.2)  and  (2.4)  imply: 


(2.6) 

(2.7) 

(2.8) 


ik  i = 2lxp» 


t = Xp. 


= ik*  •*  = ^ 


Note  that  each  A^  is  a permutation  matrix  in  view  of  (2.8)  and  the  fact  that  all 
elements  of  A^  are  0 or  1. 


Remark  2.1.  It  is  necessary  that  p < k.  From  (2.7),  since  J is  k x p,  the 
rank  of  £ is  p and  it  is  necessary  that  k 2 p.  If  k = p,  £ would  be  non-singular 
contradicted  by  (2.6).  Hence  p < k.  The  number  of  orthogonal  polynomials 
defining  the  trend  in  (2.3)  must  be  less  than  the  number  of  plots  per  block. 

3.  A necessary  and  sufficient  condition  for  a TFB  design.  The  sums  of 
squares  for  treatments  (adjusted  for  blocks  if  necessary) , blocks  (unadjusted) , 
and  trend  of  the  analysis  of  covariance  for  the  design  modelled  by  (2.5)  or 
(2.3)  are  considered.  The  general  approach  of  Searle  (1975)  is  used.  A design 
described  by  (2.5)  is  trend-free  relative  to  the  trend  effects  in  that  model 
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when  the  indicated  treatment  and  block  sums  of  squares  may  be  calculated  as  though 
the  trend  effects  were  omitted  from  the  model. 

Some  additional  notation  is  needed.  Let  the  reduction  from  the  sum  of  squares 
Y'Y  to  the  appropriate  residual  sum  of  squares  due  to  a fitting  of  the  general 
linear  model,  Y = Xy  + E,  by  least  squares  be  denoted  by  R(y).  Given  an  appro- 
priate partition  of  X and  y,  X = (X^,  and  y ' = (yj,  y^)  , the  model  becomes 
X = ♦ X2y2  + E and  R(y)  = R(Yj , y2) . Let 

(3-1)  R^X2  ^Xp  = R^Xl*  X2^  * RfYi^ 

represent  the  difference  in  reduction  of  sum  of  squares  due  to  fittings  of  the 
models,  X = &Y  + £ and  Y = X^  + E.  in  our  problem,  R(t|w,  £,  6)  and  R(jg | u , 0) 
represent  the  desired  treatment  and  block  sums  of  squares  and  R(t|p,  j3)  and 
R(£|vO  represent  the  corresponding  sums  of  squares  for  the  model  like  (2.5)  with 
trend  effects  deleted.  The  circumstances  under  which 

(3-2)  R(t|p,  £,  0)  = R(r|p,  £) 

and 

(3.3)  R(fi|w,  0)  = R(fi|y) 

are  investigated. 

The  sums  of  squares  in  (3.2)  and  (3.3)  are  obtained  by  the  Searle  method. 

We  have 

(3-4)  r(2>.  fi.  a)  - x'BjM&BjX. 

(3-5)  R(xlw,  fi)  - Y'D2XT02rD2X, 


and 
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(3.6) 


where 


(3.7) 


Rtel**.  0)  - = ^X’XgX'Y  - ^(x;y)v 


D.  = I..  - X.  (X!X. )"  X! , 
~i  ~bk  ~x» 


(3.8) 


(3.9) 


Qi  - x;d.xt,  i - 1.  2, 
h = <V  V xe),  X2  = (x,,  x0) 


and  A"  is  a generalized  inverse  of  a matrix  A,  AA  A = A.  It  is  obvious  from  (3.6) 

that  the  trend  effects  in  the  model  (2.5)  do  not  affect  calculation  of  the 

unadjusted  block  sum  of  squares  nor  do  the  values  of  the  6^.  in  . This  is 

intuitively  clear  since  trend  effects  have  zero  sums  over  plots  within  blocks 

from  (2.1)  and  the  definition  of  <#>  (t).  The  ordinary  block  sum  of  squares  cal- 

culated  from  block  totals  results  from  (3.6).  But  (3.4)  and  (3.5)  depend  on 

X and  hence  may  depend  on  the  fi*  when  the  model  (2.5)  applies.  We  show  that 
~T  JJb 

(3.5)  does  not  depend  on  choices  of  the  6*  , give  a definition  for  the  trend-free 
concept,  and  develop  a necessary  and  sufficient  condition  for  a design  to  be 
trend  free. 

Consider  Xj(X|X^)'Xj,  i = 1,  2,  in  (3.7).  Searle  (1971,  p.  20)  shows  that 
this  quantity  is  invariant  to  the  choice  of  the  generalized  inverse  and  it  may  be 
verified  that  we  may  take 


0 f Ik 
~ k ~b 


(3.10)  (XJXj)-  - J8  £1, 


0 and  (X'X2)- 
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Then 


(3.11)  and 


D.  - I..  - r XQXI  - i-  X0XI 

~i  ~bk  k ~e~B  b ~e~e 


S2  -lbIc  - 1 Me- 


Since  XJX^  = DR,  ;he  diagonal  matrix  with  i-th  diagonal  element  r^,  the  number 

of  replications  of  treatment  i,  and  X^X^  = N,  the  v x b incidence  matrix  with 

elements  n..,  n..  = 1,  0,as  treatment  i is,  is  not,  in  block  i,  Q-  and  Ql  do  not 
ij  ij  *2  *2 

depend  on  the  6*  . Since  X'Y  = T,  the  column  vector  of  trea  "^nt  totals  and 
Jt  ~x~ 

X*X  = R,  the  column  vector  with  i-th  element  r. , X'D_Y  does  not  depend  on  the 
6*  nor  does  (3.5).  Thus  (3.5)  depends  only  on  the  usual  observation  totals  and 
design  parameters  and  is  the  usual  treatment  sum  of  squares  for  the  block  design 
calculated  in  the  absence  of  trend  terms  in  the  model  (2.5). 

We  continue  with  an  explicit  definition  and  a theorem. 


Definition  3,1.  A block  design  modelled  by  (2.5)  is  trend-free  relative  to 
the  trend  in  the  model  if  R(t|u,  g,  0)  = R(t|u,  g). 

Theorem  3.1.  A necessary  and  sufficient  condition  for  a block  design  to  be 
trend-free  is  that 


(3.12)  X'X.  = 0. 

Proof.  To  prove  necessity,  equate  (3.4)  and  (3.5), substitute  from  (3.7) 
and  (3.8),  and  obtain  the  equality. 


(3.13) 


D.X  O'X'D.  = D-X  0“X'D_. 
~1~t*1~t~1  ~2~t*2~t~2 
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Pre-  and  post -multiplication  of  both  sides  of  (3.13)  by  and  X^  respectively 
reduce  the  equality  to  ^ or»  equivalently,  to 

(3.14)  X'  (D_  - D.)  X = 0. 

Use  of  (3.11)  in  (3.14)  yields  (3.12)  to  complete  the  proof  of  necessity. 

To  prove  sufficiency,  (3.12)  is  used  with  (3.11)  to  show  that  and 

that  X^D^  = and  the  equality  of  (3.4)  and  (3.5)  follows. 

Remark  3.1.  The  design  construction  problem  may  be  reformulated.  Note  that 
XJ.  = (AJ,  ...,  A^),  where  each  A..  is  a permutation  matrix.  Then  (3.12)  is 
equivalent  to 


(3.15) 


where  A^  = £a . . The  matrix  A+  has  non-negative  integer  elements  such  that  all 
row  sums  are  b and  column  sums  are  r^,  i = 1,  , v.  The  TFB  design  is  con- 
structed if  A^  satisfying  (3.15)  is  found  first  and  then  if  b permutation 

matrices  A.  are  found  such  that  Ta.  = A . 

~3  j 3 ~+ 

Notice  that  (3.15)  implies  that 


(3.16) 


b 

l 

j=l  t1=l 


I1  t 


t =1 
m 


- 0, 


for  all  i = 1,  ...,  v and  all  a e *■  . In  a TFB  design,  the  total  effect  of  each 
trend  component  over  the  plots  assigned  to  any  treatment  is  zero. 


4.  Optimality  properties  of  TFB  designs.  Let  us  suppose  that  the  primary 
purpose  of  an  experiment  is  to  compare  treatment  effects.  Then  it  is  desirable 
in  most  experiments  that  the  treatment  design  matrix  X be  chosen  such  that  each 
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possible  treatment  difference,  x^  - x^,,  i * i',  i,  i'  = 1,  v,  be  estimable 

under  model  (2.5)  as  stated  or  with  trend  components  omitted.  Designs  with 
such  a property  have  been  defined  to  be  connected.  The  selection  of  may  be 
viewed  as  a two-stage  process.  The  first  stage  is  the  determination  of  a way  of 
blocking  specified  by  the  incidence  matrix  N * X|X^  and  the  second  stage  is  the 
allocation  of  treatments  to  plots  within  blocks.  In  this  section  we  decompose 
the  class  of  connected  designs  into  subclasses  with  identical  incidence  matrices. 
Given  (2.5),  optimality  properties  possessed  by  TFB  designs  within  these  sub- 
classes are  considered. 

The  following  familiar  optimality  definitions  are  considered  over  the 
subclass  of  connected  designs  under  model  (2.5)  with  the  same  incidence  matrix: 

(i)  A design  is  V-optimal  if  the  average  variance  of  the  estimators  of  all 
elementary  treatment  contrasts,  x^  - x^, , is  minimum. 

(ii)  A design  is  A-optimal  if  the  trace  of  the  covariance  matrix  of  the 
estimators  of  any  v - 1 orthonormal  treatment  contrasts  is  minimum. 

(iii)  A design  is  D-optimal  if  the  determinant  of  the  covariance  matrix  of 
(ii)  is  minimum. 

(iv)  A design  is  E- optimal  if  the  largest  eigenvalue  of  the  covariance 
matrix  of  (ii)  is  minimum. 

4.1.  Estimability  and  connectedness.  Consider  the  model  (2.5),  first  in 
its  general  form,  and  a given  non-null  column  vector  ii  of  constants  with  proper 
dimension.  It  is  well  known  that  the  definition  of  estimability  implies  that 
u'y  is  estimable  if  and  only  if 


(4.1) 


rk(u,  X'X)  = rk(X'X) . 
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We  use  (4.1)  with  proper  selection  of  u to  determine  a necessary  and  sufficient 
condition  for  the  estimability  of  a linear  function  of  treatment  parameters  and 
consequently  for  the  connectedness  of  a design. 

The  least-squares  normal  equation  associated  with  (2.5)  is  X'Xy  = X'Y,  where 
X is  the  estimator  of  y.  Since  N'ly  = klb>  = R,  and  l^X^  = 1^  = filxp, 
it  reduces  to 


(4.2) 


where  G = X'Y,  T = X'Y,  B = X'Y,  and  W = XIY,  arc  the  observation  total  and  vectors 

~ ~r~  ~ ~0~'  ~ ~6~’ 

of  treatment,  block  and  trend  totals  respectively.  Let 


(4.3) 


B ■ (0V)<1.  iv,  - ±N,  - fex,), 


premultiply  both  sides  of  (4.2)  by  H,  and  obtain  the  reduced  normal  equation. 


for  the  estimation  of  treatment  parameters,  where 

<4-5>  £ - Er  - I®’  - kw, 

and 


(4.6) 
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The  matrix  C in  (4.4)  plays  a decisive  role  in  the  estimation  of  a linear 
function  of  the  treatment  parameters  and  in  the  connectedness  of  a design  under 
model  (2.5).  A linear  function  a'r  may  be  represented  by  u'y  with 


(4.7) 


»'  - <°»  a’.  W olxp): 


a being  a v -element  column  vector  of  constants.  It  may  be  shown  that 


(^•8)  rk (u,  X'X)  = b + p + rk(a,  C) 

by  an  argument  similar  to  that  of  Chakrabarti  (1962)  for  a model  like  (2.5) 
with  trend  effects  absent.  If  a = 0 x^, 

(4.9)  rk  (X'X)  = b + p + rkC. 

Use  of  (4.8)  and  (4.9)  in  (4.1)  demonstrates  that  a necessary  and  sufficient  con- 
dition for  the  estimability  of  a'r  is  that 

(4.10)  rk(a,  C)  = rkC. 

It  follows  further,  through  use  of  (4.10)  and  the  argument  by  Raghavarao  (1971, 
Theorem  4.2.2),  that  a design  represented  by  (2.5)  is  connected  if  and  only  if 
rkC  = v - 1.  For  a TFB  design  C reduces  to  Crt  * D_  - ^NN'.  We  summarize  our 
results: 


Theorem  4.1.  A TFB  design  under  model  (2.5)  is  connected  if  and  only  if 
• ^SS')  = v - 1.  If  a TFB  design  exists,  imposed  on  a connected  block 
design  with  a model  like  (2.5)  with  trend  effects  absent,  the  TFB  design  is 
connected  under  (2.5). 
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4.2.  Optimality  properties.  The  main  theorem  on  optimality  follows: 

Theorem  4.2.  If  a TFB  design  exists  in  a subclass  of  connected  designs 
with  a given  incidence  matrix,  the  TFB  design  is  V-,  A-,  D-,  and  E-optimal  over 
the  subclass. 

To  prove  Theorem  4.2,  some  preliminary  results,  definitions,  and  lemmas  are 
required. 

It  can  be  shown  from  (4.5)  and  (4.6)  that 

(4.11)  E(2)  = Ct  and  Var(Q)  = o2C. 

If  a'x  is  estimable,  it  is  uniquely  estimated  by  a'x  and  hence 

(4.12)  E(a't)  = a'x  and  Varfe'x)  = a2a'C+a, 

where  C+  represents  the  Moore-Penrose  inverse  of  C. 

If  any  two  square  matrices  A and  B of  the  same  order  are  considered,  we 
define  A £ B to  mean  that  A - jg  is  positive  semidefinite.  Evidently,  A i g means 
that  A is  positive  semidefinite.  Let  the  trace,  determinant  and  largest  eigen- 
value of  A be  trA,  |A|,  and  Af{ (A)  respectively.  The  following  lemmas  assist  in 
proof  of  Theorem  4.2. 

Lenina  4.1.  [Milliken  and  Akdeniz  (1977)].  If  A and  B are  two  symmetric 
matrices  of  the  same  order  and  A 2 B 2 0,  then  B*  2 A+  if  and  only  if  rk(A)  * rk(jg). 

Lenina  4.2.  For  any  two  symmetric,  positive  semidefinite  matrices  A and  B 
of  the  same  order  v and  the  same  rank  n,  if  A £ B,  then 
(i)  tr(B*)  2 tr(A+), 

(ii)  IjkbVI  2 IjcaVI, 

(iii)  Am(B+)  2 Am(A+)  and  AM(Kg+]C')  2 yKAV), 
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where  K is  any  n x v matrix  such  that  KA+K'  is  positive  definite. 

Proof.  Lemma  4.1  implies  that  B+  ^ A+  and  hence  KB+K'  * KA+K’.  The  proof 
of  (i)  is  a consequence  of  the  fact  that  each  diagonal  element  of  B+  is  not  less 

* <•  . « . •*  . I ^ iV-l'  i 

• . 1 * ■ •*  ••  ./■  • • J i . • » ' 

than  the  corresponding  one  of  A+.  The  proof  of  (ii)  follows  from  Rao  (1973,  p.  70). 
The  proof  of  (ii.)  follows  from  the  fact,  given  by  Rao  (1973,  p.  62),  that,  for 
any  symmetric  matrix  H, 

r.'Hx 
= SUP  ~Px 


where  & is  any  non-null  column  vector  of  proper  dimension. 

Lemma  4.3.  If  A is  a symmetric,  positive  semidefinite  matrix  of  order  v, 

then  n tr(M)  - 1 'Ml  > 0 for  any  n-square  principal  minor  M of  A. 

~ ~n~~n  ~ ~ 

Proof.  Since  A is  symmetric,  positive  semidefinite,  so  is  M.  Hence  M may  be 

a covariance  matrix  of  an  n x 1 random  vector,  say  £.  Let  £'  = (r.j,  ...,  ” ) 

and  M - (m„ ) . The  lemma  follows  from  the  fact  that  the  average  variance  of  all 

2 

^ i < j.  i,  j s !,•••,  n*  is  nonnegative. 

We  may  proceed  to  the  proof  of  Theorem  4.2. 


Proof  of  Theorem  4.2.  For  a connected  design,  every  elementary  treatment 
contrast,  - t.,,  i < i',  i,  i'  = 1,  ...,  v,  is  estimable.  From  (4.12),  we 
see  that 


Var (rt  - t^) 


,„ii  i ’ i * 
Cc  + c 


ii'  ♦ 

where  c , i,  i'  ■ 1,  . . . , v , is  the  (i,  i')-element  of  C . Hence,  the  average 

variance  V of  the  estimators  of  all  elementary  treatment  contrasts  is 
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_ 1 


V 7 


v_  i j,  it,  var<?i  • 

*•  -1  x<± • 


, 2 V-1  V 

(4.13)  ■ TirrT)-  I l I--11  ■ 11 

*•  AJ  i=l  i'=i+l 


♦ c 


'*•  - ii' 

4C  ) 


= ,f7^q  [v  tr(C+)  - l'C+l  ]. 
v(v  - 1)  1 ~ ~v~  ~VJ 

We  compare  the  general  matrix  C with  CQ  for  the  TFB  design  in  the  subclass 
of  block  designs  with  common  incidence  matrix  N and  an  existing  TFB  design: 

Cn  - C = k'X.XlX  ^ 0. 

~0  ~ bAir~0'o9~T  ~ 

It  follows  from  Lemma  4.1  that  C1  - C*  H and  hence 

m <mQ  IV 


(4.14) 


v tree  - Sq)  - i;(c  - Sp)^  2 0 


by  Lemma  4.3.  Applying  (4.13)  to  C and  C^,  and  using  (4.14),  we  have 


v-fr2 ! -i)  P «(£*)  - «£\l  - 7irrTTiv  tr<)  - VSjivl 


2 0, 


which  proves  V-optimality. 

Let  r be  any  (v  - 1)  * v matrix  such  that  ££•  = and  = £(v_i)xi* 

Then  [t  is  a vector  of  v - 1 orthonormal  treatment  contrasts  and  Varffjr)  * 
«Vr.  Since  £*  * Cq  and  ~£+~’  2 ££5r  > A‘»  D'»  and  E-optimality  of  the  TFB 
design  follow  from  (i),  (ii),  and  (iii)  respectively  of  Lenana  4.2. 

The  practical  effect  of  Theorem  4.2  is  to  assure  the  user  of  TFE  designs 
that  they  will  be  optimal  as  indicated  relative  to  model  (2.5)  in  comparison  with 
the  analysis  of  variance  for  the  corresponding  block  design  with  treatments  ran- 
domized over  plots  within  blocks.  The  user  of  the  TFB  designs  benefits  also  from 
the  simpler  available  analysis  of  variance  calculations. 
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5.  Analysis  of  variance  of  TFB  designs.  Calculations  and  justifications 

for  the  analysis  of  variance  of  TFB  designs  may  be  developed  by  standard  methods. 

A TFB  design  is  derived  from  a standard  block  design  for  which  the  model  is 

Y = X^y  + XgjJ  + X^t  + F,  model  (2.5)  with  trend  terms  deleted,  and  an  analysis 

of  variance  table  is  available.  That  standard  analysis  of  variance  is  used  to 

obtain  sums  of  squares  for  treatments  (adjusted  for  blocks  if  necessary)  and 

blocks  (unadjusted)  and  the  corrected  total  sum  of  squares  for  the  analysis  of 

variance  of  the  TFB  design  with  model  (2.5).  The  sum  of  squares  for  trend  with 

p degrees  of  freedom,  when  model  (2.5)  applies,  is  jV'W/b,  W = The  error 

sum  of  squares  is  obtained  by  subtraction.  The  sum  of  squares  for  trend  is 

2 

partitioned  easily  into  components,  V\L/b,  i = 1,  ...,  p,  each  with  one  degree 
of  freedom. 

While  the  paragraph  above  describes  the  necessary  calculations,  the  general 
analysis  of  variance  for  a TFB  design  is  displayed  in  Table  1.  In  that  table, 

G,  T and  jg  are  as  defined  for  (4.2)  and 


F 


±y’x’(x0x: 
k~  ~ ~B~B 


2 

The  sums  of  squares  are  independent  and,  when  divided  by  o , have  chi-square 
distributions  with  degrees  of  freedom  (d.f.)  given  in  Table  1 and  non-centrality 

parameters  obtainable  from  the  expected  mean  squares  of  the  table,  expected 

2 

mean  square  divided  by  o less  one. 

Table  1 is  simplified  in  well  known  simple  situations.  For  a TFCB  design, 
k ■ v and  r^  •«  b,  i * 1,  ...,  v,  and  then 


f - vrcib  - 
2-£^a  ■ fCC'I  - ^2>- 


So  ■ bUv  " 

I'EoftX'dv-KlpI. 


and  rk  CA  * v - 1. 
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For  a TFBIB  design  with  parameters  v,  b,  k,  r and  X, 

£o-T<l-7W'  Q'£o2  = f2'2. 

IV  * T?'«v  - PWI'  “>d  rk£o  ’ » - >• 


TABLE  1 


General  Analysis  of  Variance  for  a TFB  Design, 

Model  (2.5) 

Source  of 

d.f . 

Sum  of  Squares 

Expected  Mean 

Variation 

Square 

Blocks 

b-1 

(B'B-±G2)/k 

o2+F/(b-l) 

(unadjusted) 

D 

Treatments 

rkC_ 

02+(lr£ol)/r*£o 

(adjusted) 

Trend  Term  1 

l 

W2/b 

2 2 
o +b6j 

Trend  Term  p 

1 

W2/b 

2 2 
o *be 

P 

p 

Error 

bk-b-p-rkCg 

By  Subtraction 

2 

0 

Total 

bk-1 

Y'X-G2/bk 

— 

Remark  5.1. 

The  TFB  designs  developed  for  model  (2.5)  are  trend-free  also 

if  treatment  or  block  effects  are 

random  rather  than  fixed.  Theorem  3.1  holds 

for  random  or  mixed  effects  models 

since  condition  (3.12) 

does  not  depend  on 

distributional  assumptions  on  Y. 

Expected  mean  squares  in  Table  1 would  require 

minor  reexpression  for  a model  with  random  treatment  or  block  effects.  When 

block  effects  are 

random  and  treatment  effects  are  fixed, 

treatment  comparisons 

may  still  be  made  and  the  demonstrated  optimality  properties  of  TFB  designs 
apply. 
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6.  Examples.  Given  an  incidence  matrix  for  a block  design  and  a polynomial 
trend,  a TFB  design  may  or  may  not  exist.  Existence  theorems  and  design  con- 
struction will  be  discussed  in  a subsequent  paper.  Two  examples  are  given 
below,  one  a TF2CB  design  and  one  a TF^BIB  design.  In  the  arrays  below,  letters 
represent  treatments  and  rows  represent  blocks.  Orthogonal  trend  components, 
without  normalization,  are  given  in  the  upper  rows  of  the  arrays. 


A TF2CB  design,  v = 7,  b 


-7  -5  -3  -1  1 357 

7 1-3-5-5-317 


ABCDEFGH 

GHBACDFE 

FEHGBADC 

DCEFHGAB 


!:  A TF^BIB  design,  v = 5,  b = 10,  k = 3,  r = 6,  A = 3,  p 
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